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The Bi2Se3 class of topological insulators has recently been shown to undergo a superconduct- 
ing transition upon hole or electron doping (Cu x -Bi2Se3 with Tc=3.8°K and Pd x -Bi2Te3 with 
Tc=5°K), raising the possibilities that these are the first known "topological superconductors" or 
realizes a superconducting state that can be potentially used as Majorana platforms (L.A. Wray 
et.al, Nature Phys. 6, 855-859 (2010)). We use angle resolved photoemission spectroscopy to exam- 
ine the full details of the spin-orbital groundstates of these materials including Bi2Te3, observing 
that the spin-momentum locked topological surface states remain well defined and non- degenerate 
with respect to bulk electronic states at the Fermi level in the optimally doped superconductor and 
obtaining their experimental Fermi energies. The implications of this unconventional surface (that 
undergoes superconducting at lower temperatures) topology are discussed, and we also explore the 
possibility of realizing the same topology in superconducting variants of Bi2Te3 (with Tc ~ 5°K). 
Characteristics of the experimentally measured three dimensional bulk states are examined in detail 
for these materials with respect to the superconducting state and topological properties, showing 
that a single Majorana fermion zero mode is expected to be bound at each superconducting vortex on 
the surface. Systematic measurements also reveal intriguing renormalization and charge correlation 
instabilities of the surface-localized electronic modes. 



Three dimensional topological insulators (3D-TIs) re- 
alize a recently discovered novel state of matter which 
is not reducible to multiple copies of integer quan- 
tum Hall states (IQH), in which a topological property 
of bulk electronic states gives rise to spin-momentum 
locked two-dimensional Dirac cone surface states (a 
novel 2DEG) with a host of unusual electronic and spin 
properties^ " 1 1 1 13 i 14 . Following the axiomatic principle 
that "more is different" in many-body systems, this novel 
surface or 2DEG environment has been the subject of nu- 
merous predictions for emergent many-body states and 
device physics that cannot be readily realized by other 
material or quantum Hall-like systems^"— . Shortly fol- 
lowing the first experimental discovery of a three dimen- 
sional topological insulator, it was proposed that the in- 
troduction of superconductivity at a topological insula- 
tor (TI) surface can under the right conditions give rise 
to braidable (non-Abelian) quasiparticles which might 
be applied in quantum information science. This pro- 
posal has drawn a great deal of attention to the recent 
observation of superconducting states in doped topolog- 
ical insulators Cu^BisSes (T C =3.8°K) and Pd ;E Bi 2 Te 3 
(Tc=5.5°K)i£»2L2£. In this paper, we explore the sys- 
tematic details of the critical questions of band topol- 
ogy and electron dynamics features that constrain the 
form of superconductivity. Our results confirm our recent 
observations that the band structure of superconduct- 
ing Cuo.i2Bi2Se3 retains the well defined surface electron 



kinematics of topological insulator Bi2Se3 10 , and indi- 
cate that the surface state in superconducting variants 
of Bi2Te3 is also likely to be fully preserved at the Fermi 
level. This topologically ordered electronic system sce- 
nario creates a platform in which superconductivity is 
expected to host non-Abelian Majorana Fermion vortex 
modest - — . The observed electron kinematics are dis- 
cussed in the context of possible forms that these vortices 
may take (surface or bulk) and surprising instabilities 
revealed by Cu-doping in the topological surface band 
structure over the range x=0 to 12%. These measure- 
ments provide important clues for developing a theory 
of superconductivity in the strongly spin-orbit coupled 
electronic systems of highly topical Bi-based topological 
insulators in general. 

Early theoretical proposals for achieving non-Abelian 
surface vortices on a topological insulator followed the 
concept presented in Ref. 26 , in which superconducting 
phase coherence is introduced to a topological insulator 
surface through an interface between the topological in- 
sulator and a surface-deposited superconducting mate- 
rial. However, no superconductor has yet been found that 
can make smooth material contact without a substantial 
Schottky barrier with the known single crystal topolog- 
ical insulators such as Bi2Se3, Bi2Te3 or Sb2Te3. And 
furthermore, the presence of a superconducting material 
overlayer would in itself be an obstacle to experimental 
probing of vortices on the topological surface. There- 
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FIG. 1: (a) The crystal structure of Bi2Se3 is shown, (b) 
A low energy electron diffraction (LEED) image of supercon- 
ducting Cuo.i2Bi2Se3 taken at 200 eV electron energy shows 
a well ordered surface with no sign of superstructure modula- 
tion, (c) The hexagonal surface Brillouin zone of Cu x Bi2Se3 
is shown in red above a diagram of the three dimensional bulk 
Brillouin zone. 



fore, the very recent discovery of superconductivity with 
a significant critical temperature of Tc=3.8°K through 
Meissner effect measurements in Cu-doped Cu x Bi2Se3 2T 
was not only a material science breakthrough, but also 
popularized a new conceptual paradigm in how supercon- 
ductivity could be incorporated with topological insula- 
tor surface states for device purposes. Our measurements 
of band structure revealing the physical scenario and the- 
oretical expectations presented by superconductivity in 
Cu x Bi2Se3 are the subject of this paper. 

Undoped Bi2Se3 is a topological insulator with a large 
band gap (~300 meV) 8 , and belongs to a class of ma- 
terials M2X3 (M=Bi,Sb; X=S,Se,Te) that includes at 
least two other topological materials, Bi2Te3 and Sb2Te3, 
with smaller band gaps and more complicated bulk 
band structures^ 2 -. They share a rhombohedral crystal 
structure (see Fig. 1), with a five atom unit cell ar- 
ranged in quintuple layers and exhibit large thermoelec- 
tric power—. These materials feature almost linearly dis- 
persive spin-polarized conduction electrons on their sur- 
faces, in which electrons behave as massless relativist ic 
particles. 

Introducing superconductivity through doping of the 
TI bulk is an optimal approach to bring superconduct- 
ing phase coherence to the surface of the crystal, but 
generates a new set of challenges because the topologi- 
cal properties and crystal quality that are desired for de- 
vice applications can be directly impacted by the dopant. 
We have found that copper doping leads to increased 
variability in the crystal quality, and care is needed to 



grow and identify single-domain crystals (see LEED in 
Fig. 2(b)) with sharp quasiparticles. Of greater concern, 
numerical predictions suggest that the topological insu- 
lator surface states merge quickly with bulk electronic 
states as the chemical potential moved from the insulat- 
ing gap, and the doping required to achieve superconduc- 
tivity (10%) is far too high for well defined surface states 
to persist at the Fermi level where Cooper pairing oc- 
curs. Even though Hall effect characterization suggested 
that copper ions donated less than one carrier for each 
Cu atom 27 , extrapolated predictions based on numerics 
and band structure measurements of undoped Bi2Se3 all 
suggest that the topological insulator surface states will 
have merged with bulk electronic states after 10% copper 
doping, making the case for novel surface physics highly 
uncertain. 

However, in the case of Cua3i2Se3, nature has proven 
to be surprisingly accommodating. Direct band structure 
measurements of Cuo.i2Bi2Se3 show that the surface and 
bulk states remain separate at the Fermi level (see e.g. 
Fig. [7]), in spite of expectations to the contrary. Using 
angle resolved photoemission spectroscopy (ARPES) to 
study the material system over a range of Cu doping val- 
ues, we observe an atypical renormalization effect that 
bends the surface bands away from the bulk, and ex- 
tract band parameters that provide a framework for un- 
derstanding the novel properties that superconductivity 
may adopt in the sample surface and bulk, related to the 
spin-orbit coupled topological insulator band structure. 
Our measurements show that the conduction band of 
Cuo.i2Bi2Se3 resembles an ideal massive 3D Dirac cone, 
which is the default bulk band shape expected for topo- 
logical insulators, making the theoretical discussion and 
analytical methodology generalizable to other material 
systems. 

The ARPES measurements presented here were per- 
formed at the Advanced Light Source beamlines 10 and 
12 using 35.5-48 eV photons and Stanford Synchrotron 
Radiation Laboratory (7-22eV photons) with better than 
15 meV energy resolution and overall angular resolution 
better than 1% of the Brillouin zone (BZ). Samples were 
cleaved and measured at 15°K, in a vacuum maintained 
below 8xl0 -11 Torr. Momentum along the z axis is de- 
termined using an inner potential of 9.75 eV, fine tuned 
by the photon energy dependence shown in Fig. [8] and 
consistent with previous photoemission investigations of 
undoped Bi2Se3 8 . Large single crystals of Cu £C Bi2Se3 
were grown using methods described in Ref. 27 >. Intro- 
ducing x=0.12 copper doping for optimal superconduc- 
tivity was found to shift the z-axis lattice parameter by 
only 1.5% while leaving the in-plane lattice parameters 
and long-range crystalline order intact. Surface and bulk 
state band calculations were performed using the LAPW 
method implemented in the WIEN2K package 30 , and ac- 
curately reproduce the key elements of topological order 
in the Bi2Se3 system. Details of the calculation are iden- 
tical to those described in Ref. 8 . 

This paper is divided into three sections. In the first, 
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FIG. 2: Doping and the chemical potential: (a) The surface state of Bi 2 Se3 is a spin-textured Dirac cone. Representative 
ARPES measurements along the Y-M axis are shown for dopings spanning the range from (b) a topological insulating state to 
(f) the optimally doped superconductor. Percentage values indicate the amount by which dispersion in the upper Dirac cone 
is reduced (renormalized) relative to the undoped compound, and the Dirac point binding energy is indicated with a dashed 
white line, (g) A low resolution Fermi surface measurement over the full first Brillouin zone is shown for the superconducting 
compound, performed at high photon energy, (h) Binding energies of the surface Dirac point and bulk conduction band 
minimum are plotted as a function of doping, (i) A second derivative image of valence bands in the optimal superconductor, 
(j) Integrated ARPES intensity over the momentum range shown in (i) is plotted for different samples. 



we explore the unconventional Cu-doping behavior and 
instabilities of electrons near the surface, which are crit- 
ical to preserving the non-degenerate topological band 
structure. In the second part, we examine the topologi- 
cally ordered bulk band structure and determine param- 
eters for a material specific model of the superconduct- 
ing state. In the third part, we comprehensively analyze 
these results to evaluate how the interplay of topological 
order and superconductivity are expected to give rise to 
novel physical properties. 

I. BULK-DOPING A TOPOLOGICAL 
INSULATOR 

The effect of carrier doping on the surface of a topo- 
logical insulator is not easy to predict, as it would be 
for a generic, topologically trivial band insulator. The 
fact that surface states span a bulk band gap leads to an 
excess of conducting charge carriers on the surface, caus- 
ing charge density to be non-uniform along the z-axis. A 
non-uniform charge density leads to many-body Coulomb 
interactions that are not readily incorporated in the elec- 
tronic energies predicted by density functional theory 
(implemented here in the generalized gradient approx- 
imation, DFT-GGA). Furthermore, band energies (and 
hence carrier density) in the surface state are highly sensi- 
tive to the z-axis potential gradient at the surface because 
of the Rashba effect, adding another level of complexity. 

As a result of these and perhaps other effects, we see 
that addition of copper results in non-linear electron dop- 
ing and highly unusual changes to the Bi2Se3 surface 
state kinetics. Because of these properties, copper inter- 



calation techniques may achieve application as an exper- 
imental tool to manipulate the surface states in bismuth- 
based topological insulators such as Bi2Se3 and Bi2Te3, 
and the effects of copper doping are worthwhile to study 
in detail. In particular, the reduction in surface state 
Fermi velocity achieved by copper doping has the effect 
of increasing the energy /momentum separation between 
bulk and surface conduction bands, stabilizing the topo- 
logical surface properties as outlined later in the paper. 
Based on this, one may also speculate that the addition of 
copper or a less perturbing neutral- valence dopant could 
provide a mechanism to restore a topological surface state 
in material systems for which topological properties have 
been lost due to overlap between the bulk and surface 
bands. 

To gain a systematic understanding of the highly un- 
usual doping effect that copper has on the system, we 
present data on several stoichiometric crystals with cop- 
per added (Cu^BisSes, x=0, 0.01, 0.05, 0.12) and con- 
trast these doping levels with a case in which copper is 
substituted for bismuth as Cuo.1Bi1.9Se3. Figure Efb-f) 
shows measurements from the data sets that were used 
to calculate the surface state conduction carrier density 
(Luttinger count) for Fig. [3](d). These measurements 
were performed with high photon energies (>20 eV) and 
are used to understand qualitative material properties 
such as the approximate surface state band contour; how- 
ever we have found that lower photon energies are nec- 
essary to resolve finer details of the bulk band structure 
due to a screening effect discussed with respect to Fig. 

m 

Copper addition gradually enlarges the Fermi surface 
through electron doping and, surprisingly, brings about a 
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FIG. 3: Doping towards superconductivity leads to 
critical changes on the surface electronic and spin en- 
vironment: (a-c), Symmetrized surface state Fermi surfaces 
are displayed together with (bottom) a DFT-GGA prediction 
based on the Dirac point binding energy at x=0.05 copper 
doping. Non-zero intensity inside the panel (b) Fermi surface 
results from bulk electrons screening the surface state, as de- 
scribed in Fig. 6. (d) The total number of charge carriers in 
the bulk and at the surface is calculated from the Luttinger 
count ( bz "area > x ^ ^ OT tne doubly degenerate bulk band). 
DFT-GGA predictions show the carrier density obtained by 
aligning DFT band structure with the experimentally deter- 
mined binding energy of the Dirac point. 



strong reduction in surface electron velocities through re- 
shaping of the surface state conduction band. The slope 
of the surface state band represents particle velocity, and 
is reduced by an increasing amount as the Cu concentra- 
tion is raised to 5%, by which point electron velocities 
near the Dirac point are 30% slower than at the Dirac 
point of undoped Bi2Se3. The surface electron veloc- 
ities recover slightly at superconducting doping but are 
still approximately 20% slower than in the undoped com- 
pound. Higher binding energy valence band structure 
evolves through a rigid shift as doping is changed (Fig. 
[2jj)), suggesting that with the exception of adding elec- 
trons, copper does not significantly alter the bulk elec- 
tronic properties or perturb the more localized (tightly 
bound) electronic states. The energy spacing of peak 
features at binding energies greater than 1 eV in the 
integrated density of states are essentially independent 
of doping, implying that renormalization of the surface 
state conduction band is not due to a change in the bulk 
material Hamiltonian such as the strength of spin orbit 
coupling. 

As the surface state Fermi surface enlarges due to 
added carriers it also becomes hexagonally anisotropic, 
as outlined in Fig. [3ja-c). Such hexagonal-like deforma- 
tion makes the surface susceptible to spin-fluctuation or 
magnetic instabilities^. Tunability of the surface state 
kinetics and anisotropy is important for control of un- 
conventional ordered states that may appear uniquely in 



topologically ordered materials. Carrier density in the 
surface state is much greater than density in the bulk, 
as estimated by the Luttinger count (Fig. Efd)), suggest- 
ing that the surface may carry a screened negative charge 
within the normal state of the superconductor. Attempt- 
ing to force the creation of Cu#i replacement defects by 
adding less bismuth results in very weak hole doping for 
Cuo.1Bi1.9Se3 (Fig. |2pi)), raising the bulk conduction 
band entirely above the Fermi level so that the mate- 
rial is a traditional topological insulator. These results 
establish that copper atoms can add holes or electrons 
depending on their net preference for occupying different 
kinds of site in the Bi2Se3 lattice. Copper atoms interca- 
lated between van der Waals-like bonded selenium planes 
(Ciiint) in the crystal are thought to be single electron 
donors, while substitutional defects in which copper re- 
places bismuth in the lattice (Cu Bi) contribute two holes 
to the system^. The bulk doping levels observed in our 
data suggest that the relative number of copper atoms 
occupying these two configurations evolves from a ratio 



of 



Cm. 



3.5 at x=0.01 to just 2.2 in superconducting 



x=0.12 crystals. 

Reducing the slope of the surface state increases sep- 
aration between the bulk and surface conduction bands, 
a characteristic that preserves stability of the topolog- 
ical order. In Fig. Efb), the undoped Bi2Se3 surface 
dispersion is overlaid on the band structure of x=0.12 
0^612863, demonstrating that with no reduction in the 
undoped band velocity, it is not expected that the surface 
and bulk bands will remain separate. The calculation in 
Fig. HJc) suggests one possible mechanism by which the 
surface band velocity could be reduced. Increasing the 
spacing between the outermost two atomic layers by just 
0.2 A from the bulk equilibrium yields a much "pointier" 
surface Dirac cone that is in better agreement with the 
experimentally observed shape, and increases the upper 
Dirac cone dispersion by 16% while lowering the energy 
of states in the lower Dirac cone. One explanation for 
reduced velocity may therefore be that the crystal sur- 
face of undoped Bi2Se3 is more strongly relaxed (lower- 
ing the valence band energy), but as more electrons are 
added to the conduction band it becomes energetically 
favorable to lower the conduction band energies by re- 
ducing the interatomic spacing to a value closer to its 
bulk equilibrium position. Other possible explanations 
and contributing factors exist and will require extensive 
experimental characterization as mechanisms for tuning 
the surface properties of topological insulators. 

The surface of a topological insulator is a complex en- 
vironment for photoemission measurements due to the 
vertical-axis anisotropy and strong topological protection 
of the Dirac-like surface states. Although DFT-based cal- 
culations predict only a single surface conduction state, 
surface sensitive ARPES measurements performed at 
high photon energy show an additional distortion local- 
ized above the bulk conduction band in energy, which 
appears to occur separately from the topological sur- 
face conduction state. When the band is measured with 
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FIG. 4: Doping towards superconductivity leads to 
changes on the surface band (and Fermi) veloci- 
ties: (a) Photoemission measurements at high photon energy 
(E>20 eV) are shown for non-superconducting Cu x Bi2Se3. 
Particle velocities in the surface state upper Dirac cone are 
reduced by 30% after the addition of x=0.05 copper doping. 
The energy scale of panel (a,right) is offset by a rigid shift to 
align the Dirac point binding energy with the panel at left, (b) 
The surface state dispersion from undoped Bi2Se3 is overlaid 
in green on a photoemission image of Cu cc Bi2Se3 at super- 
conducting composition (x=0.12), measured with low energy 
(9.75 eV) photons to observe (blue) the bulk dispersion. By 
linear extrapolation, we project that if the surface state ve- 
locity were not reduced after Cu doping, the bulk and surface 
bands would intersect below the Fermi level at the energy 
labeled with a dashed white line, (c) When a 12-quintuple- 
layer slab calculation is modified by shifting the outermost 
(Se) atomic layer 0.2 A away from the next (Bi) atomic layer 
(layer spacing increased from 1.641.84 A), dispersion in the 
upper Dirac cone increases by 16%. This calculation demon- 
strates that the surface state dispersion is very sensitive to 
minor changes in the surface chemical environment. 



higher bulk sensitivity (Fig. [5fc,left)), it can be clearly 
resolved as a narrow positive-mass paraboloid, however 
the inside of the paraboloid becomes filled in with inten- 
sity when higher photon energies are used (h^>20 eV, 
see Fig. [5jc,right)). Figure [2] contains further examples 
of filled-in bulk band structure. 

A likely explanation for this phenomenon relates to the 
many-body Coulombic interaction between electrons at 
the surface and in the bulk, which is not accounted for 
by DFT. As we have observed through the systematic 
doping measurements summarized in Fig. [3fd), the sur- 
face state carrier density of Cu^Bi2Se3 and Cuo.1Bi1.9Se3 
is much greater than the carrier density in the bulk, and 
holds a negative charge. Theory suggests that surface 
state penetration into the crystal is approximately given 
by a ratio of the surface state band velocity and the bulk 
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FIG. 5: Bulk electrons screen the surface charge: (a) 

The crystal structure of CuccE^Ses, viewed parallel to the 
(111) cleavage plane. Copper can donate electrons when in- 
tercalated between Se layers, or may accept electrons if it 
substitutes for Bi. The bulk and surface state electron den- 
sity distributions are sketched as discussed in the text, (b) 
The large surface state carrier density is screened by the 
bulk band, causing "bulk" states to be hole doped (occur 
at smaller binding energy) in close proximity to the surface, 
(c) ARPES measurements with (left) large and (right) sub- 
nanometer penetration depth are compared. 



band gap (v^/A), which is less than 1 nm for Cu^E^Sea 
(3eVA/0.4 eV = 7.5A for x=0.12). Given this scenario, 
the bulk and surface carrier densities (and charge den- 
sity) are expected to conform to a distribution similar 
to that drawn in Fig. [5fa). We note that details of the 
absolute charge density distribution may depend on the 
broader sample environment and boundary conditions. 

Panel EJb) shows an approximation of hole doped 
states the bulk electrons may occupy to screen the surface 
state within 1 nm of the surface. These states constitute 
a continuum of upwardly shifted "band images" that fill 
in the the inside of the bulk dispersion due to (poten- 
tially quite complex) Coulomb interaction with the sur- 
face state. Because the screening effect is highly local- 
ized at the surface, these hole doped band images have 
broken translational symmetry along the z-axis, and are 
generally somewhat visible at arbitrary k z values (e.g. 
see Fermi surface maps in Fig. [3fa-c)). In a classical 
picture based on the Poisson equation, it is expected 
that screening will occur over less than 0.5 nm 32 , but 
this is clearly a non-physical model for detailed predic- 
tions on such a short length scale because of the discrete 
nature of the lattice, derealization of the surface state 
and numerous related concerns. Qualitatively, reduction 
in the screening length as carrier density increases is ex- 
pected to strengthen hybridization between bulk and sur- 
face electronic states, and may be partly responsible for 
the reduced surface state band velocities we observe. 

Topological insulator surfaces are often unstable over 
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FIG. 6: Time evolution of Cu-substituted 

Cuo.1Bi1.9Se3: Unlike Cu-intercalated Cu^Bi^Ses The 
chemical potential in Cu-substituted Cuo.1Bi1.9Se3 slowly 
rises over time after the surface is exposed to vacuum. The 
binding energy of the surface state Dirac point is indicated 
with an arrow for spectra measured (a) 5 minutes (b) 30 
minutes (c) 1.5 hours and (d) 9 hours after cleavage. The 
undoped Dirac point binding energy is indicated with a 
dashed white line in panel (d). Alignment of crystalographic 
axes in (a) differs slightly from other panels, (e) The charge 
density distribution expected near the surface of a highly 
insulating, freshly cleaved Bi2Se3 sample is shown, (f) Two 
mechanisms that can potentially cause the surface chemical 
potential to shift upwards over time are the migration of 
positively charged selenium vacancy defects attracted by 
the negatively charged surface state, and the adsorption of 
positive ions from the vacuum (mostly H + ). 



a long period in vacuum, becoming increasingly elec- 
tron doped over a period of minutes or hours (see e.g. 
RefjLI 2 .). Of the Cu-doped samples studied, only the 
surface of Cuo.1Bi1.9Se3 was observed to change over a 
24 hour time period, with a gradually increasing chemical 
potential as shown in Fig. [6ja-d). The chemical poten- 
tial rises quickly in the first hour and a half, with ap- 
proximately 0.002 electrons added per surface Brillouin 
zone, however it takes another 7.5 hours to add 0.002 
more electrons, suggesting that the effect saturates or 
will only occur for highly insulating samples. Theories 
suggest that this time dependence is unrelated to the fun- 
damental topological order of topological insulator band 
structure 3 -, but it nonetheless constitutes a noteworthy 
behavior of the known topological insulator compounds. 
Recent experimental studies in Ref. 11 have shown that 
strongly perturbing the surface energetics of a topologi- 
cal insulator results in the formation of new surface bands 
and Dirac points, in configurations that respect the topo- 
logical insulator order. 

Analyses have shown that the rising surface chemi- 
cal potential is probably due to gradual accumulation 
on the surface of positively ionizing gas molecules 33 , and 
it is well established that depositing positive ions on the 



FIG. 7: Band topology at superconducting composi- 
tion: (a) Momentum dependence of the bulk and surface 
conduction bands in superconducting Cuo.i2Bi2Se3 is mea- 
sured through the 3D Brillouin zone center with low energy 
(9.75eV) photons for enhanced bulk sensitivity, (b) Disper- 
sion along the Y-M and Y-K directions is traced from the 
blown-up measurements shown in panels (c-d), revealing sig- 
nificant anisotropy in the surface state dispersion. 

crystal surface will donate negative charge carrier o 7 i n i 12 . 
The data presented here suggest one additional instabil- 
ity that may lead to a similar effect, namely that the 
negative charge we have observed in the surface state 
will attract positively charged selenium vacancy defects, 
causing the surface to gradually become more strongly 
electron doped as summarized in Fig. EJe-f). Such an 
effect would weaken greatly as the surface chemical po- 
tential approaches the bulk band. The electric field from 
the surface is screened over a distance similar to the ratio 
of surface and bulk charge densities, which may be tens 
of nanometers or longer for highly insulating samples 32 . 

These systematic doping studies show an unexpect- 
edly large spectrum of surface electronic behavior, from 
hexagonal warping of the Fermi surface to strong renor- 
malization of band velocities, all of which occurs while 
the chemical potential changes by an unexpectedly small 
amount. These features highlight the strong sensitivity 
of the surface state kinetics to intrinsic properties such as 
the relaxation of lattice spacing near the cleaved crystal 
surface and complicated interactions between the surface 
and bulk electrons. 



II. BULK ELECTRON KINETICS AND 
SUPERCONDUCTIVITY 

The compounds Cu^Bi2Se3 and Pd x Bi2Te3 are the 
first representatives of a novel class of superconducting 
materials with strong spin-orbit coupling. Unlike con- 
ventional superconductors or noncentrosymmetric super- 
conductors with spin-split bands, these materials have an 
inversion symmetric crystal structure with doubly degen- 
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FIG. 8: Band topology along the Y-M axis: (a) DFT- 

GGA calculations predict that the surface and bulk conduc- 
tion band structure will overlap first (at highest binding en- 
ergy) along the Y-M axis. An inset shows the predicted over- 
lapping Fermi surface at the experimentally derived chemical 
potential, with the k z -projected bulk band shaded in blue. 

(b) The experimentally determined Fermi surface contours 
are shown, with black dots marking specific Fermi momenta 
identified from panel (d) and Fig. 7. (c) Dispersion along 
the z-axis is examined by varying the incident photon energy. 

(c) Dispersion along Y-M is shown for momenta spanning the 
bulk k z dispersion to establish that there is no degeneracy 
between the bulk and surface state electrons at k z ^0. 



erate bulk bands. In this case, spin-orbit coupling leads 
to strong interband mixing that cannot be captured by 
any single-band model, and is the origin of the topologi- 
cal insulator state. As a starting point for studying bulk 
superconductivity in Cu^E^Sea, the spin-orbit coupled 
character of its band structure needs to be examined in 
detail. 

Doping into the superconducting regime (x=0.12) 
raises the Fermi energy to 0.25 eV above the bulk conduc- 
tion band minimum, placing the Fermi level in a highly 
linear regime of bulk band dispersion (Fig. [7J). There- 
fore, Cooper pairing in the superconducting state is char- 
acterized by the properties of these linearly dispersive 
electronic states close to the chemical potential. Bulk 
band structure with a "relativistic" lineshape that tran- 
sitions from parabolic to linear as momentum increases 
is a direct consequence of band inversion required for 
the topological insulating state, and has also been ex- 
perimentally observed in topological insulator Bii-^Sb^ 
alloys 5 -. In Bi2Se3 band inversion occurs at the T-point 
(k x =k y =k z =0) in the center of the bulk conduction band 
leading to the expectation of a spin-orbit induced Dirac- 
like bulk band. In our experiments this is indeed the 
case - the bottom of the conduction band for electrons in 
the bulk is found at the three dimensional T-point, which 
also lies inside of the upper surface state (SS) Dirac cone 



without touching it. Bulk Fermi momenta of 0.110±3 
A -1 and 0.106±3 A -1 are observed along the Y-M and 
Y-K directions respectively. Carefully tracing the band 
(Fig. EJb-d)) reveals a Fermi velocity of 3.5 eVxA along 
Y-M and 4.1 eVx A along Y-K, estimated within 50 meV 
of the Fermi level. The gap between bulk valence and 
conduction bands appears to be unchanged upon copper 
doping, consistent with the observation from Section I 
that the spin orbit coupling strength and bulk Hamil- 
tonian appear to be unchanged. These experimentally 
determined bulk and surface Fermi surfaces are convex, 
while the numerically predicted Fermi surfaces are con- 
cave (Fig. [8(a-b)), emphasizing the critical importance 
of experimental investigations to understand the band 
structure and electronic states of topological insulators. 
The numerics can be approximately reconciled with ex- 
perimental measurements by multiplying ("renormaliz- 
ing") DFT-derived band energies by a factor of 2 above 
the bulk conduction band minimum. 

The kinetic behavior of electrons with momentum per- 
pendicular to the cleaved surface is examined in Fig. 
[8](b-d) by varying incident photon energy. The z-axis 
Fermi momentum observed in Fig. [8(c) is 0.12±1 A -1 , 
suggesting that the bulk electron kinetics are three di- 
mensionally isotropic with only slightly reduced veloci- 
ties in the out-of-plane direction. The essential kinetics 
of a positive-mass band with less than 1% carrier oc- 
cupancy as is the case here can generally be evaluated 
from measurements that pass through the band mini- 
mum along the three high-symmetry directions (k^, k y 
and k z ). However, band structure calculations suggest 
that bulk dispersion along the Y-M axis will intersect 
the surface state band before any other location (Fig. 
[Ha)), and may be less-than-parabolic for k 2 ^0, prompt- 
ing a closer investigation. Band structure measurements 
shown in Fig. [8(d) examine dispersion along the Y-M 
axis at out-of-plane momenta spanning the bulk k z dis- 
persion, confirming that the separation between bulk and 
surface bands increases monotonically as momentum is 
increased along the z-axis. Exploring the full momentum 
space reveals that the surface and bulk electronic states 
at the Fermi level are well separated by 0.04 A -1 in mo- 
mentum and an energy spacing of about A#=130 meV, 
with closest proximity in the k z =0 plane. 

A minimal two-band model for the band structure of 
Cu^Bi2Se3, which correctly takes into account of spin- 
orbit coupling and crystal symmetries, has been proposed 
recently 1 —. The k • p Hamiltonian near Y in the 3D Bril- 
louin zone is given by 

H (k) = mcf x + v(k x a z s y - k y a z s x ) + v z k z a y , (1) 

Here a matrices are Pauli matrices acting on a two-band 
basis. The eigenvalues of <j x = 1(— 1) label the con- 
duction (valence) band at k = 0, which is mainly an 
even (odd) linear superposition of two p z orbitals on the 
top and bottom Se-Bi bilayer within the five-layer unit 
cell. The parameter s represents Pauli matrices labeling 
electron spin, and k z is along the [111] direction per- 
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FIG. 9: Massive 3D Dirac fermions in Cii^E^Ses: 

(a) Energy- momentum dispersion of bulk electrons through 
the 3D BZ center is compared with Dirac-like (v=6 eV-A) 
and classical (parabolic) fits with the same effective mass 
of 0.155M e . Surface-state (SS) dispersion is plotted with a 
massless linear fit. (b, blue) The bulk band dispersion ap- 
proximated from photoemission data is plotted at k 2 =0, sur- 
rounded by (green) a massless "light-like" Dirac dispersion 
with identical Dirac velocity (v c ), illustrating that bulk elec- 
tron kinetics are in the relativistic Dirac regime, (c) A small 
panel shows the difference between measured bulk band en- 
ergy and the Dirac-like fit curve in panel (a). 

pendicular to the cleavage plane of the rhombohedral 
crystal. As shown in (pQ), the spin-orbit coupling term 
is off-diagonal in band basis. The form of Ho resem- 
bles the Dirac Hamiltonian for massive electron and po- 
sition in 3 + 1-dimensional QED. This becomes transpar- 
ent when we introduce a set of 4 x 4 Gamma matrices: 
T cr x , Ti a z s y , T 2 ~cr z s x , T 3 a yi so that 
Ho = raTo + v(k x Ti + k y T2) + v z k z Ts- We note that the 
velocity of the massive Dirac fermion is anisotropic due 
to lack of rotational symmetry in the crystal. 

We now fit the band dispersion measured by ARPES 
with that of the model Hamiltonian Ho: 

E± = ±^m 2 + v 2 kl + v 2 k 2 + v 2 z k 2 z . (2) 

The curves in Fig. [9fa) demonstrate that the conduc- 
tion band can be closely described using a Dirac mass of 
m=0.155 M e (electron masses) and Dirac velocity of v=6 
eV-A. Contrasting this fit with a classical (parabolic) 
dispersion of identical mass reveals that electron kinet- 
ics begin to enter the linear relativistic regime within 



~0.1 eV of the Fermi level Fig. [9ja-b). A slight bend 
in the dispersion centered near 90 meV binding energy 
(Fig. Efc)) may suggest strong electron-boson interac- 
tions in the system also ubiquitously observed in other 
superconducting materials 34 or a shift in the balance of 
orbital and spin-orbit terms defining the electron kinet- 
ics. This model effectively addresses symmetry mixing 
from spin-orbit coupling and the qualitative shape of 
the conduction and valence bands, but does not capture 
fine details of the valence band dispersion. Experimen- 
tal and numerical studies have shown that the bulk va- 
lence band of topological insulators in the M2X3 family 
is warped 8,12 , with a local minimum (rather than maxi- 
mum) at the Brillouin zone center due to orbital-kinetic 
terms dominating over spin orbit interactions to reverse 
the sign of electron velocity for a small range of momenta. 
Spin-orbit symmetry mixing captured by Equations 1-2 
is very weak in the immediate vicinity of the T-point, al- 
lowing for significant deviations of this sort, but quickly 
becomes the dominant kinetic term at larger momenta 
where Cooper pairing takes place. 

The observed spin-orbit band structure kinematics de- 
fine several key parameters that give insight into aspects 
of the superconducting wavefunction without the need 
for a detailed model. The average Fermi velocity ob- 
served from the measurements shown in Fig. [Tfb-d) 
suggests that the superconducting coherence length is 
approximately 2000A (f ~ 0.2 x hv F /K B T c = 0.2 x 
3.8eVA/(K B x 3.8°K) = 2000A), assuming minimal 
scattering and a superconducting gap related to Tc by a 
BCS-like perturbative Cooper pairing mechanism. This 
value is about 1-2 orders of magnitude greater than the 
correlation lengths estimated from band structure of the 
cuprates (£ ~ 100-200A) or cobaltates (£ ~ 200A)2^. 
Correlation length and carrier mobility set the phase or- 
dering temperature scale, which can thus be estimated 
to be T e - 60,000°K (T e = £0 x h 2 n e /2m*= 2000A x 
ft 2 x 10 20 cra- 3 /(2 x .155Me) = 60,000°K), four orders of 
magnitude larger than the superconducting critical tem- 
perature of 3.8°K. These properties of large coherence 
length and high phase ordering temperature suggest that 
the superconducting phase transition and bulk supercon- 
ducting gap formation can be understood within a mean 
field picture, although to describe the spin-orbit symme- 
tries and collective properties of the topologically ordered 
superconducting electrons one needs a material specific 
model as we will discuss further in Section III. Based 
on mean field considerations, we expect the supercon- 
ducting gap to be about ~0.6 meV (3.5xK jB Tc/2=0.6 
meV), much smaller than the state-of-the-art resolution 
of ARPES. 

The spin textured topological surface state is key to 
many of the most interesting properties that can emerge 
from the interplay of topological order and broken sym- 
metry. Therefore for any superconducting doped topo- 
logical insulator it is critical to experimentally identify 
if after doping the surface state is fully preserved at the 
Fermi level where Cooper pairing occurs, as our mea- 
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( a ) Bi2Te3 GGA ( b ) p-type Bi 2 Te 3 (c) n-type Bi 2 Se 3 




FIG. 10: Band topology in superconducting composi- 
tion of Bi2Te3: (a) The band structure of Bi2Te3 is pre- 
dicted by the generalized gradient approximation (GGA) for 
a six quintuple- layer slab, with lighter colors marking states 
with greater occupancy of the outermost quintuple layer, (b- 
c) The band structures of Bi2Te3 and Bi2Se3 measured by 
ARPES with high photon energies (>20eV) are shown side 
by side, and an inset displays the Bi2Te3 Fermi surface in a 
±0.2A _ 1 window, revealing only the surface state. Panels 
at the bottom show band structure at higher binding energy, 
with dashed lines highlighting a pattern of surface resonance 
that is common to topological insulators in the M2X3 series. 
The spectrum in panel (b,top) was measured immediately af- 
ter cleavage, before any significant observed shift in the sur- 
face chemical potential. 



surements shown in Fig. [7] and Fig. [8] have established 
to be the case for Cuo.i2Bi2Se3. As identified earlier in 
this Section, the separation of bulk and surface electronic 
states is 0.04 A -1 in momentum and A#=130 meV in 
energy, defining the scale on which perturbations can be 
applied for device engineering without compromising the 
capability to form collective electronic states that are lo- 
calized entirely at the surface. For the surface state to be 
well defined, it must be fully separate (non-degenerate) 
from the k^-projected bulk band continuum at the Fermi 
level. A scenarios in which this condition is not met is 
demonstrated by the numerically predicted Fermi surface 
in Fig. [8fa), and the important role played by unusual 
surface instabilities to achieve a non-degenerate surface 
state in the 0113312863 material matrix has been dis- 
cussed in Section I. 

Now that superconducting states with large transi- 
tion temperatures have been identified in chemically 
(Pd ;E Bi2Te3, Tc=5.5 °K) 28 and mechanically strained 
(Tc~3 °K) 36 versions of the closely related compound 
Bi2Te3, it is desirable to perform the same sort of band 
structure analyses for that system. Numerical simula- 
tions such as that shown in Fig. [TOTa) suggest that 
the surface state upper Dirac cone at the superconduct- 
ing chemical potential will be well defined and found at 
smaller momenta than the bulk conduction states. How- 



ever, photoemission measurements (Fig. HDT b)) have yet 
to reveal any sign of the bulk conduction bands at the 
Fermi level, in spite of a significant p-type Hall carrier 
concentration of lxl0 19 cm -3 . The carrier density of p- 
type as-grown Bi2Te3, shown in Fig. HDTb), is thought 
to be similar to superconducting variants of Pda3i2Te3, 
which cannot yet be grown with a high degree of crys- 
talline homogeneity. A region of photoemission intensity 
from 0.1 to 0.3 eV binding energy merges with the lower 
Dirac cone, and is blurred in a pattern that resembles the 
surface-screened bulk conduction band of Bi2Se3, but has 
no visible intersection with the Fermi level. 

As with undoped Bi2Se3, the measured band structure 
of Bi2Te3 is characterized by larger group velocities than 
are predicted by first principles calculations. The mea- 
sured T-M Fermi velocity of the surface state is 1.9 eVA, 
as compared to 2.3 eVA in the slab calculation. Surface 
resonance features at larger binding energies such as the 
diamond shape traced in the bottom panels of Fig. [TOTb- 
c) are shortened on the energy axis by a factor of ~2 in 
the Fig. [TUTa) calculation, similar to the observed renor- 
malization of the Cua3i2Se3 conduction band. With no 
definite experimental observation of the bulk conduction 
states, and significant quantitative uncertainty in first 
principles calculations, it is currently challenging to de- 
termine whether or not the bulk and surface states of 
superconducting Bi2Te3 overlap. 

The bulk band structure of topological insulators has 
previously been studied to only a limited degree because 
of the screening effect described in Fig. [5] and the bulk 
insulating character, which limits experimental investiga- 
tion of the conduction band. Due to the low photon ener- 
gies used and electron doped chemical potential, the re- 
sults presented in this section give a clear view of the bulk 
conduction band kinetics of Bi2Se3, a compound that has 
been termed a "hydrogen atom" topological insulator for 
its simplest-case realization of topological insulator 4,8 . 
These band parameters, including a nearly isotropic mas- 
sive Dirac dispersion and the complete separation of sur- 
face and bulk conduction states provide the key details 
of electron kinetics and topology to formulate a theory 
of the novel superconducting properties of Cuo.i2Bi2Se3. 
As we will discuss in Section III, superconductivity at the 
surface of Cuo.i2Bi2Se3 cannot be conventional given the 
observed band properties, and is expected to be the first 
realization of one of two novel forms of superconductivity 
characterized by electronic topology. 



III. IMPLICATIONS: NON-ABELIAN 
SUPERCONDUCTOR OR TOPOLOGICAL 
SUPERCONDUCTOR? 

Our experimental observations have several impor- 
tant implications for new superconducting states in 
Cu^Bi2Se3. First, our data show that the spin-textured 
surface and doubly degenerate bulk bands are well sepa- 
rated. This establishes that bulk superconducting pair- 
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U/V 



FIG. 11: Surface Majorana vortices vs. topological 
superconductor The "m//i" ratio derived from ARPES pa- 
rameters is critical to defining whether the ground state will 
be a (A2-state) topological superconductor with odd parity 
superconductivity, or (Ai-state) will host non-Abelian Majo- 
rana Fermion surface vortices with an even parity supercon- 
ducting wavefunction. Parameters "U" and "V" represent 
matrix elements for intraorbital and interorbital scattering 
amplitudes from the superconducting pairing interaction, as 
described in the text and Ref.— . 



ing occurs in the presence of a non-degenerate, spin po- 
larized two dimensional topological surface state. As a 
result, the proximity effect with the bulk superconduc- 
tor is sufficient to induce the pairing of surface states, 
without the need of using an external superconductor 
as proposed earlier 26 . The resulting 2D superconduct- 
ing state on the surface is highly nontrivial: it is fully 
gapped but supports zero-energy Majorana bound states 
localized in the vortex core. These Majorana states have 
non-Abelian statistics and are potential building blocks 
of a topological quantum computer. For this reason, we 
call the superconducting state on the surface "a 2D non- 
Abelian superconductor". Moreover, since Cu^E^Sea is 
a type-II superconducto r 2 ^ 28 , the required vortex on the 
surface naturally appears at the end of a vortex line in 
the bul k 37 i 38 , and can be easily generated by applying a 
magnetic field (Fig. fT3H)). 

Theoretical studies suggest that the number of Majo- 
rana fermions at the end of a vortex line is partly de- 
rived from bulk band structure properties at the Fermi 
leve l 38 1 39 . If the Fermi level is inside or very close to the 
bulk band gap of a topological insulator, then an odd 
number of Majoran fermions (e.g. one) is expected be 
bound to each end of a vortex line. If the Fermi en- 
ergy is raised far above the conduction band minimum 
of the topologically inverted band gap, the ID vortex 
line will eventually cross a topological phase boundary 
beyond which the end of each vortex will bind an even 
number of surface Majorana modes (e.g. zero or two). 
A recent study based on DFT band structures has pre- 
dicted that the transition will occur when the Fermi level 
is raised higher than 0.24 eV above the conduction band 
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FIG. 12: Superconductivity on a topological surface: 

(a) A phase diagram compares (ARPES expt) the mea- 
sured superconducting topology with preliminary expecta- 
tions based on cases in which (A) each Cu atom donates one 
doped electron, (B) the experimental chemical potential is 
applied to GGA band structure, and (C) the experimental 
chemical potential is applied to the band structure of un- 
doped Bi2Se3. (b) States within 5meV of the Fermi level are 
shaded for the k 2 =0 plane of (center) a fully gapped doped 
topological insulator and (right) one example of a topological 
superconductor. A panel at left shows a scenario in which de- 
generacy between the bulk and surface electrons allows super- 
conductivity to be achieved without novel topological proper- 
ties. 



minimum 39 ; however, the factor of ~2 discrepancy we 
have observed between calculated DFT band kinetics and 
experimentally measured band structures in Cu x Bi2Se3 
and Bi2Te3 implies that the transition occurs at a higher 
energy, possibly Em~2x0.24 eV ~0.5 eV above the con- 
duction band minimum in Cu :E Bi2Se3. Because we ob- 
serve the Fermi level to be only ~0.25 eV above the 
conduction band minimum for optimally doped super- 
conducting Cuo.i2Bi2Se3, the expectation based on our 
measurements of bulk band structure is that only one 
Majorana mode will be found at the end of each vortex. 
The well denned surface state observed in the measure- 
ments presented here is expected to give rise to a single 
Majorana mode that is localized within several nanome- 
ters of the crystal surface^, and our ARPES measure- 
ments clearly rule out the possibility of a second Majo- 
rana bound state arising from surface states. If a second 
mode exists (contrary to present predictions), it must 
come from bulk electrons, and therefore will be extended 
far more deeply along the vortex line into the material. 
In either case, the observation of non-Abelian phenom- 
ena made possible by the Majorana modes will depend 
on physical parameters that have yet to be investigated, 
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such as detailed interactions and tunneling effects be- 
tween vortices. 

In addition to these surface properties, our data also 
provide important clues for understanding superconduc- 
tivity in the bulk of Cu^E^Sea. Strong spin-orbit cou- 
pling in the form of 3D massive Dirac bands plays a key 
role in the pairing symmetry of Cu^E^Ses. As Ref^ 
shows, an unconventional spin- triplet pairing whose order 
parameter is odd under spatial inversion strongly com- 
petes with the conventional spin-singlet pairing. Specif- 
ically, two phenomenological parameters U and V were 
introduced in the two-orbital model, which denote the 
intra- and inter-orbital electron interactions respectively. 
It was found that triplet pairing wins at small electron 
density, provided that inter-orbital attraction exceeds 
over the intra-orbital one (see phase diagram in Fig. 
[TT]) . In reality, this condition is reasonable because larger 
intra-orbital repulsion tends to significantly weaken the 
phonon-mediated attraction. 

The key band structure parameter constraining the 
form of superconductivity in the perturbative model pre- 
sented in Section II is the ratio of rest mass (band cur- 
vature) to the doped chemical potential (rn/fi). In an 
ideal Dirac band structure described by the Hamiltonian 
in Eq. (1-2), the rest mass is equal to half of the band 
gap between the conduction and valence bands, or 150 
meV. Because the observed chemical potential is 250 meV 
above the conduction band minimum, this gives a ratio 
of — ~ I (blue arrow in Fig. [TT]) . which is small enough 
to allow the rare theoretical possibility that the bulk su- 
perconducting wavefunction will have an odd (-1) parity 
symmetry value. However, the warped shape of the va- 
lence band (partical-hole asymmetry) that we have dis- 
cussed above causes different mass terms to be needed 
to describe the conduction and valence band lineshapes, 
with a larger mass term needed for the conduction band. 
For the experimentally derived fit parameters plotted in 
Fig. [9j the ratio is raised to — ~ §, which would suppress 
the likelihood of odd superconducting pairing. The pres- 
ence of a "kink" in the dispersion centered near 90 meV 
binding energy causes a more linear dispersion locally 
at the Fermi level (see Fig. [9fc)), making the velocity of 
Cooper-paired electrons consistent with the lighter band- 
gap rest mass and suggesting that the ratio of ^ ~ ^ 
best describes the superconducting electron kinetics. The 
"mass gap" between the valence and conduction bands 
comes directly from spin orbit coupling interactions, and 
is therefore a direct indicator of how orbital symmetries 
are influenced by spin-orbit coupling. Different detailed 
models may shift the boundary in the Fig. [TT] phase dia- 
gram, but are not expected to change the overall physical 
scenario in which both even and odd states are possible 
for differently tuned pairing interaction strengths. 

As a result, triplet pairing wins for a large and reason- 
able range of interaction parameters. Interestingly, the 
triplet superconducting state realizes a novel time rever- 
sal invariant "topological superconductor" phase, which 
is the electronic analog of the BalianWerthamer (BW) 



phase in superfluid He-3. Recent theoretical research 
has shown that three dimensional topological insulators 
doped into a superconducting phase are classified by an 
integer topological invariant ("n"), just as the undoped 
topological insulator state is identified by a Z2 invari- 
ant (parity invariant)^. The topological superconduc- 
tor state expected under odd-parity pairing in Cu x Bi2Se3 
is identified by a non-zero value of "n". A hallmark of 
this topological phase is gapless surface Andreev bound 
states within the pairing gap (see Fig. [T2T b. right) and 
Fig- fT3t c)). These Andreev bound states are unrelated 
to the electron surface states of the parent topological in- 
sulator Bi2Se3. The bound states host Bogoliubov quasi- 
particles which are two-dimensional Major ana fermions 
propagating on the surface. However, these delocalized 
Majorana modes are unlike the Major ana Fermions local- 
ized at vortices in the topological surface state, and no 
proposals currently exist for their use in quantum com- 
puting. 

While triplet superconductivity has been postulated 
in several material systems (e.g. Sr 2 Ru04, UGe2, ZrZn 2 
and URhGe 40 45 ), it is associated with the spontaneous 
breaking of time reversal symmetry in all of these cases 
due to the magnetic moment of triplet Cooper pairs. 
Breaking of time reversal symmetry rules out any three 
dimensional topological superconductor stated. The 
spin-orbit coupled band structure of Cu x Bi2Se3 allows 
a novel form of triplet Cooper pairing in which the total 
spin and orbital angular momentum is (see Ref^ and 
the inset of Fig. [TT]) . unlike the known p + ip pairing 
in which Cooper pairs have orbital angular momentum 
corresponding to a magnetic moment. The odd parity 
superconducting state suggested by the band parameters 
reported here for Cuo.i2Bi2Se3 is therefore non-magnetic, 
and the odd parity "A2" superconducting state in the 
Fig. [12] phase diagram is identified as a topological su- 
perconductor. 

In summary, our measurements show that surface 
instabilities cause the spin-helical topological insulator 
band structure of Bi2Se3 to remain well defined and non- 
degenerate with bulk electronic states at the Fermi level 
of optimally doped superconducting Cuo.i2Bi2Se3, and 
that this is also likely to be the case for superconduct- 
ing variants of p-type Bi2Te3. These surface states pro- 
vide a highly unusual physical setting in which super- 
conductivity cannot take a conventional form, and is ex- 
pected to realize one of two novel states that have not 
been identified elsewhere in nature. If superconducting 
pairing has even parity, as is nearly universal among the 
known superconducting materials, the surface electrons 
will achieve a 2D non-Abelian superconductor state with 
non-commutative Majorana fermion vortices that can po- 
tentially be manipulated to store quantum information. 
Surface vortices will be found at the end of bulk vortex 
lines as drawn in Fig. [T3H). If superconducting pairing 
is odd, the resulting state is a novel state of matter known 
as a "topological superconductor" with Bogoliubov sur- 
face quasi-particles present below the superconducting 
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FIG. 13: A Majorana platform: (a) Topologically protected surface states cross the Fermi level before merging with the 
bulk valence and conduction bands in a lightly doped topological insulator, (b) If the superconducting wavefunction has even 
parity, the surface states will be gapped by the proximity effect, and vortices on the crystal surface will host braidable Majorana 
fermions. (c) If superconducting parity is odd, the material will be a so-called "topological superconductor", and new states 
will appear below Tc to span the bulk superconducting gap. (d) Majorana fermion surface vortices are found at the end of 
bulk vortex lines and could be manipulated for quantum computation if superconducting pairing is even (Ai). 



critical temperature of 3.8 °K. As drawn in Fig. fl3t c). 
these low temperature surface states would be gapless, 
likely making it impossible to adiabatically manipulate 
surface vortices for quantum computation. The unique 
physics and applications of the topological superconduc- 
tor state are distinct from any known material system, 
and will be an exciting vista for theoretical and experi- 
mental exploration if they are achieved for the first time 
in Cua;Bi2Se3. 

Comprehensive spectroscopic images measured over a 
wide range of photon energies reveal the detailed massive 
Dirac-like bulk and surface conduction states and their 
topological structure at optimal superconducting com- 
position, and reveal that the critical electronic param- 
eter ratio ^ is sufficiently large to render either parity 
state of superconductivity possible depending on the bal- 
ance of Cooper pairing interactions present. These elec- 
tronic and spin-orbital parameters provide the essential 
framework for device development based on the unique 
topological surface physics, and are an essential guide 



for future exploration of superconductivity in the strogly 
coupled spin-orbit systems of topological insulators. 
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